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' Abstract: We present an overview of information theory approach (both in its extensive and nonextensive versions) applied to high 

energy multiparticle production processes. It will be illustrated by analysis of single particle distributions measured in proton-proton, 
>^ proton-antiproton and nuclear collisions. We shall demonstrate the particular role played by the nonextensivity parameter q in such 

analysis as summarizing our knowledge on the fluctuations existed in hadronizing system. 
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1 Introduction 

Information theory has been estabhshed to be extremely useful tool in many 
branches of science |[T] . Its first application to multiparticle production processes is 
very old [2], however, it was used only sporadically since then PH]. Here we shall 
present its methods and analyse, as example, single particle spectra obtained in 
multiparticle production processes at high energy [Hlllj . Such processes are respon- 
sible for the majority of observed cross section. If /ci 2 are four-momenta of colliding 
objects and W = ^/s = \J {ki + /c2)^ is invariant energy of collision, then the mean 
multiplicity of produced secondaries increases with energy like {N) oc In^W or 
(A^) oc W"' (with a ~ 0.5) reaching in heavy ion collisions values of the order of 
~ 10^-^ ~ 10^. This fact makes analysis of such processes very difficult, in par- 
ticular it is impossible to describe them from first principles. On the other hand, 
understanding these processes is very important because they provide background 
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for more sophisticated measurements and they are themselves valuable source of 
information on the hadronization process in which initial energy W converts into 
observed secondaries (mostly mesons tt and K). One is looking therefore for phe- 
nomenological models and first candidates were statistical models of different kinds 
[5,6j. With necessary modifications and improvements they are actually in com- 
mon use also nowadays and are very successful [7j. They all assume some kind 
of local thermal equilibrium forming at early stage of collision process. However, 
this assumption is not necessary because the characteristic "thermal-like" (i.e., 
exponential) behaviour of some observables arises always whenether one restricts 
itself to observation of only a small part of some large system |Bj . This is precisely 
situation encountered in multiparticle production processes where usually only a 
small part of produced secondaries is registered by detectors and out of them only 
very limited number is subjected to final analysis. The necessary averaging over 
all unmeasured degrees of freedom introduces therefore a kind of randomization 
or a heat bath, action of which can be summarily described by a single parame- 
ter r, a kind of " effective temperature" of the usual thermodynamics . It could, 
however, happen that such "thermal bath" is more complicated (exhibiting, for ex- 
ample, some intrinsic fluctuations or collective flows) [9]. In such cases parameter 
T must be supplemented by additional parameter summarily describing the action 
of these new factors - in this way the nonextensivity parameter q appears (such 
"heat baths" are not extensive anymore) |2j and notion of nonextensive statistics 
enters in a natural way |[TD] . 

Why use information theory methods in analysing multiparticle production pro- 
cesses? To illustrate this let us suppose that experiment measures some new dis- 
tribution of secondaries. Immediately theoreticians rush to describe the new data 
and provide a number of models, which differ drastically in what concerns their 
assumptions, nevertheless they all do explain experimental results in a remarkable 
good way. What has happened? To answer this one has to introduce notion of infor- 
mation: data considered contain only limited amount of information and all models 
possess it as well. They differ therefore in what concerns additional redundant in- 
formation they posses, which reflects not so much the state of the measured object 
but rather the state of minds of scientists proposing these models. To quantify 
this one has to resort to methods of information theory and introduce informa- 
tion entropy as a measure of information PI21I31- The problem one is facing can 
be formulated in following way. Suppose one has system S, depending on some 
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variable x G {xmin, Xjnax), on which one performs finite number n measurements 
providing us with Ri^,,,^n results. The question is: how to obtain the most plausible 
{least biased) and model independent probability distribution, p{x), of variable x? 
The answer is : look for such p{x) (normalized to unity) which maximizes in- 
formation entropy under constraints given by results of measurements. To account 
for many known features of the system under consideration (like long-range corre- 
lations, fractal-like structure of the phase space or intrinsic fiuctuations, which will 
be of special interest to us later \T1\) we shall use Tsallis entropy jTDj characterised 
by parameter which for q ^ 1 becomes the usual Shannon entropy [T] 
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under constraints^ 

"-^TnCLX 

dxPS'''>(x) ■ [p(x)]'' = (flC^'), = 



(2) 



As result one gets 



P{x) = ^expg 



(3) 



where 



X' 

A. 



X 



1/(1-9) 



=^ exp(-). 



(4) 



Here are Lagrange multipliers obtained by solving constraint equations and Zq 
ensures normalization of p{x) to unity. It should be mentioned here that fiuctua- 
tions in the parameter 1/ A of the exponential distribution lead to the g-exponential 
with q given by normalized variation of the parameter 1/A piTUT^ . 

* For our limited purposes this form is adequate. Formalism using escort probability distributions, Pi = 
Pi/J2iPi': leads to distributions of the type c[l — (1 — q)x/l]'^^^^~'^\ which is formally identical with @, 
c[l - (1 - g)x/L]^/(^"'3), provided we identify: Q = l + {q-l)/q, L = l/q and c = {2-Q)/L = l/l. Now 
(x) = L/(3 - 2Q) = 1/(2 - q) and < Q < 1.5 (to be compared with 0.5 < g < 2). Both distributions 
are identical and the problem, which of them better describes data is at this level of sophistication not 
important. 
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Long time ago it was shown using this method (with q = 1) \T\ that in or- 
der to fit single particle distributions dN/dy in rapidity variable y = |ln(E' + 
Pl) / l'n.(E — pl) and multiparticle distributions P{N; N) measured at that time 
it was enough to assume that: {i) transverse momenta of produced secondaries 
Pt are limited (i.e., phase space is effectively one dimensional) and {ii) not the 
whole available energy ^/s is used for production of the observed secondaries but 
only its fraction K G (0, 1) called inelasticity, the rest of energy is taken away by 
the so called "leading particles". Closer scrutiny of numerous models competing 
at that time revealed that all of them contained these two assumptions (explicitly 
or implicitly) and that these assumptions consisted the only part common to all 
of them. All other additional assumptions were spurious and unnecessary and as 
such could be safely abandoned without spoiling agreement with data. 



2 Results 



We shall present now examples of analysis of single particle distributions in ra- 
pidity for pp and pp collisions at energies varying between ^/s = 20 GeV to 1800 
GeV jlj and first results for central nuclear collisions (see Fig. [T]). Our input con- 
sists of: {i) the energy of collision, Ecm = y/s, out of which, according to result of 
[2], we shall retain only fraction Kq for production of secondaries; {ii) the mean 
multiplicity of charged secondaries produced in nonsingle diffractive reactions at 
given energy: flch = —7.0 + 7.2s'^-^^^ (out of which we construct the total number 
of produced particles, = ^fich needed in our procedure); (ii) the fact that phase 
space is essentially one dimensional with mean transverse momentum slightly en- 
ergy dependent given by (pr) = 0.3 + 0.044 In (-^5/20) (all secondaries will be 
assumed to be pions of mass fi = 0.14: GeV) and {iv) that multiplicity distribu- 
tions of produced secondaries is not Poissonian (like it was in [21) but of Negative 
Binomial (NB) type 



Here E is energy of the observed particle produced in collision process and is its longitudinal 
momentum, i.e., projection of its total momentum on the collision axis. The respective transverse momenta 
Pt are integrated over and enter only via their mean value, {pt), defining the so called transverse mass 
of the particle: fix = + (pt)'^ with fi denoting its real mass. Notice that E = //rcoshy whereas 
PL = HTsinhy. 
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where parameter k {k > 1) affects width of P{n) and is given by its normahzed 



-0.104 + 0.058 In ys. (6) 



variance 

1 cr^ (rich 



k {uchY (rich) 



This last point is new and most important because it can be accounted for only in 
nonextensive version of our method with q > 1. This is because the value of k^^ 
may be understood as the measure of fluctuations of mean multiplicity caused by 
some dynamical factors. When there are only statistical fluctuations in hadronizing 
system then P{n) is of Poissonian type. Intrinsic (dynamical) fluctuations would 
mean fluctuating mean multiplicity n. In the case when such fluctuations are given 
by a gamma distribution with normalized variance D{n) then, as a result, one 
obtains the NB multiplicity distribution with 1/k = D{n) (because convolution of 
Poissonian and Gamma distributions results in NB one |T7j). Assuming now that 
these fluctuations contribute to nonextensivity, i.e., that they define parameter 
g = 1 + D{n), one is lead to the use of g > 1 (i.e., nonextensive formalism) with 
g = 1 + where k is parameter defining NB form of P(n), energy dependence 
of which is given by. eq. 



With such input the most probable rapidity distribution is 

]5(2/) = = ^exp^(-/3q • /ircosh?/) , (7) 



where Zq and Pq are not free parameters, like in similar formulas used in applica- 
tions of statistical models but are given by the normalization condition and 
by the energy conservation constraint (yUy = \/Jj^~+~{pt)^ , see ji] for details): 

j dyfiT ■ cosh?/ • [p{y)Y = \ ■ (8) 

-Y 



Symmetry of the collision (identical particles) implies that momentum is conserved 
in this case automatically. There are therefore two parameters, nonextensivity q 
responsible for fluctuations (as mentioned above) and g'-inelasticity Kq to be de- 
duced directly from data essentially in a model independent way, cf. Fig. [TK. There 
is, however, question concerning physical meaning of K,q for g > 1 case considered 
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here (notice that our q changes from 1.05 to 1.33 when going from = 20 GeV 
to 1800 GeV whereas the "partition temperature" Tq = l/Pq changes, accordingly, 
from 1.76 GeV to 62.57 GeV The point is, as was stressed in [19j, that in 
case of q > 1 our energy Kq contains also interaction between particles and is not 
directly seen in rapidity distribution. However, the quantity (inelasticity) 




1 

y p, [GeV] y 



Fig. 1. Examples of applying eq. (jSJ- (o) fit to rapidity spectra for charged pions produced in pp and 
pp collisions at difi'erent energies |14| . (6) fit to pT spectra from UAl experiment ^^I- (c) fits to rapidity 
spectra for negatively charged pions produced in PbPh collisions at different energies (for the most central 
— 7% interval of impact parameter of collision jl6j). The best fit for 17.3 GeV is for two-source case, cf. 
Fig. |2f, with (7 = 1 and Ay = 0.83. Total inelasticity in this case is K = 0.58. See text for other details. 

is the correct fraction of energy used for production of secondaries particles (it 
happens to be fairly constant with ^/s varying slightly around Kq ~ 0.5 |jlj). 
Fig. [IId shows example of similar fit but this time to the observed pr spectra at 
energies 200, 500 and 900 GeV [15j (using eq. (Ej) with energy iircoshy replaced 
by transverse momentum px)- As can be noticed fits are very good in the whole 
range of pr with the corresponding values of Tt = 1/Pt (in GeV) and qr equal to 
{TT;qT) = (0.134; 1.095), (0.135; 1.105) and (0.14; 1.11) for energies 200, 500 and 
900, respectively. These values must be compared with the corresponding values 
of (T = Ti^q = qi) for rapidity distributions, which are equal to, respectively: 
(11.74; 1.2), (20.39; 1.26) and (30.79; 1.29) at comparable energies gj. The question 
then arises, what is the actual value of the parameter q describing fluctuations of the 



6 



mean multiplicity, i.e., affecting, according to the previous discussion, multiplicity 



distributions? It is sensitive to p = \lp\ + and, as we have seen here, both pi 
and pt show traces of fluctuations by leading to q > 1 . At present we can offer only 
some approximate answer because there are no data measuring pt distributions at 
all values of rapidity i.e., providing correlations between parameters {T;q) = 
{Tl'i Ql) for longitudinal momenta (rapidity) distributions and (T; q) = [Tt] qr) for 
transverse momenta distributions. Noticing that q — l= cr^(r)/r^ (i-e., it is given 
by fluctuations of total temperature T) and assuming that cr^(r) = a'^{Ti)+a'^(TT), 
one can estimate that resulting values of q should not be too different from 

^ = r2 — + 1' (10) 



which, for Ti ^ Ty, as is in our case, leads to the result that q ~ q^, i.e., it is 
given by the longitudinal (rapidity) distributions only. 

Finally, on Fig. [Tb we show fits to the recent NA49 data IE] on tt^ production in 
PbPb collisions at three different energies (per nucleon). The values of {q; Kg = 
3 • i^'f")) obtained are: (1.2; 0.33) for 17.3 GeV, (1.164; 0.3) for 12.3 GeV and 
(1.04; 0.22) for 8.6 GeV. The origin of q > 1 in this case is not yet clear. The 
inelasticity seems to grow with energy. It is also obvious that for higher energies 
some new mechanism starts to operate because we cannot obtain in this case 
agreement with data using only the energy conservation constraint. The possibility 
offered in this matter by some closer inspection of properties of eqs. (El) and ((Tj) are 
displayed in Fig. [21 Whereas upper panels are self-explanatory, the lower demand 
some attention. They show how p{y) can be distorted and it could be that one of 
the mechanisms presented there is showing up in high energy NA49 data mentioned 
above. Three mechanisms are listed there (only q = I case is considered). Fig. 2d 
shows the changes in p{y) introduced by ^/-dependent px, especially when {px) 
grows substantially towards y = {in such a way as to keep the averaged {px) the 
same as original 0.4 GeV/c). In this case one observes increase of p{y) at small 
values of y but it seems to be too wide to explain the NA49 data. 

Fig. [2h shows result of the special kind of momentum dependent residual interac- 
tions discussed in [20]. In our case it would result in 

p{y) = i • exp [-/3 ■ iiT coshy - X ■ iit\ sinhy|] (11) 
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Fig. 2. Examples of properties of eqs. © (VF = 200 GeV, N = 60, {pt) = 0.4 GeV/c). (a) spectra obtained 
for different parameter q practically coincide with extensive spectra produced by masses M* = M/{3 — 2q). 
For all remaining panels q = 1. (6) quasiscaling for spectra obtained for composition of smaller masses 
producing respectively smaller number of secondaries with M/N =const. (c) spectra with different values 
of M/N ratio differ drastically, (d) the pT growing towards the center of rapidity phase space has dramatic 
effect on spectra, (e) the role of momentum dependent residual interaction (cf. eq. (jllj) ) is visualised. (/) 
the effect of two sources separated in rapidity by 2Ay with combined energies and masses equal M is 
visualised. 

with additional constraint imposed now on the modulus of the momenta: 



J dy iiT\si^^y\p{y) = 1^ ■ — 

Y 



M2 



N 



(12) 



Here W = K ■ ^/s and M = J2i yfW^I^i with pu = fiTsinhyi (sum is over all 
produced particles, i = 1,...,A^, and k, is parameter telling us the fraction of 
momentum turned into interaction). As it can be seen in Fig. Eh, the effect of this 
interaction although visible is very weak and probably negligible at present data. 

The most promising seems to be the conjecture that invariant energy W = K • ^/s 
hadronizing into particles consists in reality with two subsources of energies 
Wi = which are displaced from y = by some rapidity Ay (being a free 
parameter here) and producing N/2 particles each, Wi^2 = W/[2cosh.{Ay)]. As one 
can see, the obtained shape is very similar to what is observed experimentally (cf. 
Fig.Hb). 
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3 Summary 



In some types of data on multiparticle production processes one frequently encoun- 
ters ambiguity concerning the question, which of the particular models used at that 
time is the correct one [2j . Such situation arises always when data contain only lim- 
ited amount of information. To select this information one has apply information 
theory methods, which are widely known and used in other branches of science [T|. 
The examples shown here show that information theory ideas can be successfully 
used also to analyse data from multiparticle production processes and that in this 
way one gets highly model independent estimation of some quantities, in our ex- 
ample it was inelasticity parameter K [4j. Actually, in [21] we have attempted to 
fit single particle distribution without a priori introducing neither inelasticity nor 
fluctuations in mean multiplicity and found that it is possible only with q < 1. 
The reason turned out to be simple: in this case the most important factor was 
decreasing of the available phase space to mimic the action of the inelasticity and 
this can be done only with q < 1. The fit was not as good as shown here with notion 
of inelasticity introduced explicitly but it was not very bad either. As in pTUT^ 
we are stressing here the connection between necessity to use nonextensive version 
of information theory and some intrinsic fluctuations existing in the hadronizing 
system. Finally, as was clearly displayed in Fig|2l our method seems to be also 
very useful in describing the gross features of the single particle spectra observed 
in heavy ion collisions. In particular, it seems that with growing energy of collision 
there is room for some new mechanisms, not present at more elementary nucleonic 
collisions. This point deserves therefore some special scrutiny in the future. 

We would like to finish mentioning that there exists also example of very successful 
use of the information theory approach to describe not only single but also double 
particle spectra, namely the so called Bose-Einstein correlations observed between 
identical bosons in all multiparticle production data, see [22j. 
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